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We present a precise determination of the polarizability and other proton structure dependent contri- 
butions to the hydrogen hyperfine splitting, based heavily on the most recent published data on proton 
spin dependent structure functions from the EG1 experiment at the Jefferson Laboratory. As a result, the 
total calculated hyperfine splitting now has a standard deviation slightly under 1 part-per-million, and is 
about 1 standard deviation away from the measured value. We also present results for muonic hydrogen 
hyperfine splitting, taking care to ensure the compatibility of the recoil and polarizability terms. 



I. INTRODUCTION 

In this article, we consider precision calculation of the 
hyperfine splitting (hfs) of hydrogen, with the goals of cal- 
culating the hfs to a part-per-million (ppm) accuracy for 
ordinary (electronic) hydrogen and of extending the calcu- 
lation to the muonic hydrogen case. 

Experimentally, the hfs of the hydrogen ground state is 
known to 13 significant figures in frequency units [l|, 

E h fs(e~p) = 1 420.405 751 766 7(9) MHz . (1) 

On the theoretical side, at the level of a ppm accuracy, 
the QED corrections are not in question. Rather, achieving 
the stated accuracy requires better evaluation of corrections 
from the finite size of the proton. Finite size corrections 
come from two-photon exchange, Fig. [1] where there is the 
possibility that the photons are individually hard and can 
see deeply into the proton. For one photon exchange, pro- 
ton structure plays no role at the ppm level because the 
momentum transfer is necessarily very low. 

Presently, our ability to numerically deal with quantum 
chromodynamics (QCD), the theory that governs how mat- 
ter is bound together to form a proton, is insufficient to cal- 
culate proton structure corrections to the desired accuracy. 
Instead, the corrections can be related to proton structure 
information measured in electron-proton scattering. The 
information is codified in terms of Pauli and Dirac form fac- 
tors F\{Q 2 ) and i 7 ^*? 2 ) for the elastic case and structure 
functions g\(v, Q 2 ) and giiy-, Q 2 ) for the spin-dependent in- 
elastic case. Here Q 2 = —q 2 , where q is the 4-momentum 
transferred from the electron and v is the energy trans- 
ferred from the electron in the lab frame; one can also use 
x — Q 2 /(2m p v), where m p is the proton mass. 

Recently reported @j data on gi{v,Q 2 ) from JLab spur 
the present study. The data are from the EG1 collaboration, 
and extend the measurements of giiv.Q 2 ) down to much 
lower Q 2 than previously available. The relations between 
the hfs and g\(v, Q 2 ) weight heavily on the low Q 2 data, so 
the latest data, which include a careful analsys of statistical 
and systematic errors, lead to a more accurate and reliable 
hfs calculation. 

The proton structure dependent corrections can be di- 
vided into Zemach, recoil, and polarizability corrections, to 
be defined shortly. The first two depend entirely on elas- 
tic intermediate states in the two-photon exchange, and all 



contributions from inelastic intermediate states are in the 
third. Our main, though not exclusive, focus will be upon 
the polarizability corrections, which contain the dependence 
upon gi(v,Q 2 ) and which have had larger statistical and 
systematic uncertainties limits than the other two terms. 
The situation has now improved and we will find that the 
uncertainty in the polarizability corrections is now compa- 
rable to the uncertainty in the elastic structure dependent 
terms. Hence we shall evaluate these also, using up-to-date 
form factors, and discuss the uncertainty limits in the cal- 
culations of all terms. 

For muonic hydrogen hfs, there is currently no measure- 
ment, but one may be possible [1], so it is appropriate to 
quote a calculated result. In the muon case, the struc- 
ture dependent corrections are in total larger than the QED 
corrections, because the former have a lepton mass pro- 
portionality, while the latter are to a first approximation 
independent of the lepton mass. 

The muon case prompts a discussion of the definitions 
of the recoil and polarizability corrections. The sum of all 
proton structure corrections is unambiguous. However, the 
separation between the recoil and polarizability corrections 
depends upon a protocol. The issue is that the elastic and 
inelastic corrections separately have (after an overall lepton 
mass, me, is factored out) logarithmic divergences in the 
nig — > limit. For convenience, the polarizability correc- 
tions have been defined Q by taking the inelastic correc- 
tions and adding an elastic-looking term to cancel the loga- 
rithmic singularity. An identical term is subtracted from the 
recoil corrections, and the overall sum is unchanged. 

The term added to form the polarizability correction must 
satisfy the criteria that it cancel the existing me, — > di- 
vergence, and that it introduce no new divergence. This 
does not uniquely fix the residual non-divergent part of the 
term. For the electron case, the choice is standard. Hence, 
one can in principle add calculations of electronic hydrogen 
polarizability and recoil corrections from different sources 
without worry. For the massive lepton case, it appears that 
there are two different protocols, which agree in the me — > 
limit but not otherwise. Hence, there is a need for care 
in combining muonic hydrogen calculations from different 
sources, or else for a unified calculation of all the proton 
structure dependent terms, as we do here. 

Our calculations and results are detailed in Section QU 
The relevant formulas are first summarized and discussed, 



followed by numerical evaluations for the electronic and 
muonic hydrogen systems. Section [TTT| summarizes our con- 
clusions. 



II. FORMULAS AND CALCULATIONS 
A. Formulas and calculations 

The calculated hyperfine splitting can be given as @, @] 

E his (t-p) = (1 + A QED + A£ vp 

+ A£ vp + A^+A s )££, (2) 
where lepton H~ is either e~ or [i~ and the Fermi energy is 
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Mass m r = mgm p /(m p + mt) is the reduced mass and 
is the Rydberg constant (in frequency units). By conven- 
tion, the Bohr magneton \ib is inserted for the lepton and 
the measured magnetic moment fi p is used for the proton. 
The constants on the right-hand-side are well enough known 
to evaluate the Fermi energy to 0.01 ppm. 

The first four corrections are due to QED, hadronic vac- 
uum polarization, muonic vacuum polarization, and weak 
interactions (Z° exchange). They are well enough known 
not to require discussion here. 

The proton structure dependent corrections are 



A s = A z + Al + A 



pol 
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The subscripts stand for "Zemach," "recoil," and "polar- 
izability." The measured value of -Ehfs(e~p) and calcu- 
lated values of other quantities implies a "target value" 
As = — 32.77 ± 0.01 ppm for ordinary hydrogen [1, @. 

The structure dependent corrections can all obtained by a 
dispersive calculation of the two-photon exchange diagram 
(Fig. [1]), as pioneered by Iddings [7] and by Drell and Sul- 
livan [4j. The reason for separating the result into three 
terms is partly to shorten individual formulas and partly for 
historical reasons. We only quote the results, reserving a dis- 
cussion of the derivation, particularly for the massive lepton 
case, for a later report. 

We start with the polarizability corrections. They are 
usually given in the limit m e —> @, 0, H H [HJ- To our 
knowledge, the only previous exception is in the article of 
Cherednikova, Faustov, and Martynenko Including the 
lepton mass, our result for the polarizability corrections is 
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FIG. 1: Upper part: the full box; lower part: the box with 
elastic intermediate states only. 
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where n p is the proton anomalous magnetic moment in nu- 
clear magnetons, x t h = Q 2 / '{2m p m. n + m 2 + Q 2 ) with m w 
the charged pion mass, and 
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The plain /?x,2 auxiliary functions, introduced by De 
Rafael [9], are 



#i(t) = -3t + 2t 2 + 2(2-t)Vt(t + 1), 
&(t) = 1 + 2r - 2v/r(r + 1) , (8) 
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Finally, F 2 {Q 2 ) is the (elastic) Pauli form factor, normalized 
byF 2 (0) = Kp . 

The polarizability terms come mainly from inelastic in- 
termediate states in Fig. [TJ hence the appearance of the 
structure functions gi 2 . The term containing F 2 is the term 
described in the introduction, which is inserted to cancel the 
divergence that appears in the gi term of Ai in the massless 
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lepton limit. As mi — > 0, one can show that 71.2 — > /?i.2 
and P\(ti) — > 1, and further, for Q 2 — > 0, /3 x (t) — * 1. Then 
the Gerasimov-Drell-Hearn [12, Ej| sum rule, 

lim — f / dx 9l (x,Q 2 )=-K 2 p (11) 

both ensures that the second term of Ai diverges in the 
massless limit, and that the first term will regularize it. 

Our polarizability correction agrees with [ll] for the gi^ 
terms, which are unique. Some choice is possible for the 
F 2 terms, and we and [ll[ made different choices. Further 
explanation of this point joins the discussion of the recoil 
correction, below. 

The Zemach corrections are fl4| 

A z = ~2am r r z (l + S r z ad ) , (12) 
where rz is the Zemach radius 

1 



The electric and magnetic Sachs form factors are 

G M (Q 2 ) = F 1 (Q 2 )+F 2 (Q 2 ) 1 

G E (Q 2 ) = F 1 {Q 2 ) - ^F 2 (Q 2 ), (14) 

and the Dirac form factor is normalized with -Fi(O) = 1. 

The extra radiative correction <5^ d is given in [X5L fl6|] . For 
the dipole form factor, Ge.m(Q 2 ) oc (1 + Q 2 /A 2 )~ 2 , one 
finds 5% d = {a/37:) [2 ln(A 2 /m 2 ) - 4111/420] = 0.0153, 
using the standard value A 2 = 0.71 GeV 2 . For other form 
factors that we use, the changes in 6 r ^ d have a 0.01 ppm 
or smaller effect upon the hyperfine splitting. 

The leading order recoil corrections are (l5l. 
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Factoring out an overall me, there remain recoil terms 
that diverge like ln(m^) as mi — > 0; hence the mg — > limit 
is not taken. Further, and in contrast to the Zemach cor- 
rections, the recoil corrections are not zero in the static and 
pointlike proton limits. [The static limit neglects the Q 2 de- 
pendence of the form factors, so that Fi{Q 2 ) — > F\(0) = 1 
and F2(Q 2 ) — > ^(0) = n p ; the pointlike limit additionally 
takes k p — > 0.] Thus, part of the recoil correction is struc- 
ture independent. However, that they are overall structure 
dependent is clear, and so it is proper to include them here 
along with Az and A po i. 

Notice that the last term in the recoil correction is the 
negative of the F 2 term from the Ai polarizability correc- 
tion. These are the terms that were added and subtracted 
to ensure that Ai contained no divergence in the massless 
lepton limit. 

We specify the term here using a historical criterion. An 
alternative non-dispersive calculation of the elastic contribu- 



tions alone, the lower part of Fig.[TJ inserts photon-proton- 
proton vertices 

r„ =l v F 1 {Q 2 ) + ^-a up qPF 2 {Q 2 ), (16) 

for incoming q, and does the loop integral directly. In mod- 
ern times, one should hesitate to do the calculation this way 
because there is no reason to think the vertex representation 
is correct when the intermediate proton is off shell. (The 
dispersive calculation is not subject to the same criticism be- 
cause it obtains the real part of the two-photon corrections 
from the imaginary part, which only requires knowing the 
vertices when all protons are on-shell.) However, the direct 
calculation is historically older than the dispersive one, and 
is still often quoted; for relatively modern source see (l7j . 
It is possible to choose the F 2 terms in the polarizability 
and recoil corrections to cancel the zero mass divergence in 
one case and give the historical result in the other, and that 
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is the choice we made. Ref. which showed only the 
polariza bility term, differs from us in the F| terms in A po i 
and this can be traced to a different choice early on [18]. 



B. Ordinary hydrogen polarizability corrections 



For electronic hydrogen, take the van 
the well-known result @, 0, H H, 0, 
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Information on gi(x,Q 2 ) is obtained from polarized elec- 
tron on polarized proton inelastic scattering, and the low- 
est Q 2 data come from the EG1 experiment at JLab, re- 
ported in 0, [3 with data details posted on the High 
Energy Physics database at Durham University (UK). The 
EG1 data have beam energies of 1.6 and 5.7 GeV, and give 
gi{v,Q 2 )/Fi{v,Q 2 ) [Fi^ 2 {v, Q 2 ) are the spin-independent 
structure functions] for 28 Q 2 bins with central values from 
0.0496 to 4.96 GeV 2 , and with W ranging from threshold to 
about 1.65 GeV for the lower -Ebeam data and from thresh- 
old to about 3.1 GeV for the higher -Ebeam data. We ob- 
tain Fi(v,Q 2 ) in the resonance region from the Christy and 
Bosted parameterization [20], and in the scaling region from 
combining the F 2 (v, Q 2 ) fit of the NMC collaboration [2l| 
with the fit to R (the ratio of longitudinal and transverse 
photon cross sections) from SLAC E143 [22j]. Where there 
is no EG1 data, we use fits to g\ from Simula et ai [23| 
in the resonance region and from SLAC E155 [24] in the 
scaling region. 

For Q 2 below Q\ = 0.0452 GeV 2 (the lower edge of the 
lowest bin) there is no data, and we complete the integrals 
by interpolating data between higher Q 2 and zero Q 2 . For 
B\, which is proportional to Q 2 as Q 2 — > 0, this is possible 
because the GDH sum rule [H, 13 f' xes tne slope, so that 
for small Q 2 , 
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We obtain c 1B = 4.94 ± 0.30 ± 1.22 GeV~ 4 by fitting to 
B 1 (Q 2 ) in the data region below Q 2 = 0.3 GeV 2 . (This is 
somewhat larger than the c x we quote in [25], partly because 
we are here fitting Bi instead of I\, but more because of 
improvements in the data.) The contribution to Ai from 



TABLE I: Contributions to A po i for electron case. Statistical, 
systematic, and modeling errors, in that order, are given in 
the parentheses and discussed in the text. 
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the low Q 2 range is thus 
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where rp is the Pauli radius of the proton, coming from the 
expansion of F 2 (Q 2 ). 

The structure function g 2 gives a small contribution to 
the hfs, because the auxiliary function (3 2 is small for the 
kinematics where the g 2 integral has its main support. This 
is fortunate, since g 2 for the proton is not well measured. 
The A 2 's in our Tables are based on a model for g 2 provided 
by the EG1 collaboration, which we also used in [25j and 
which is heavily based on the MAID parameterization [26| 
of existing photo- and electroproduction data. Given the 
overall lack of data that is specific to g 2 , we assigned 100% 
error limits to the A2 determnation; even so, the contribu- 
tion to the overall uncertainty of A po i is not large. A more 
detailed discussin of A2 is given in Sec. Ill El 

Results for the polarizability correction are broken down 
in Table [T] for one particular parameterization of the elastic 
form factor F 2 (Q 2 ). In this Table, "systematic errors" mean 
systematic errors that come from the listed data @, [3, and 
"modeling errors" come from error limits accompanying the 
parameterizations [23] that we use to complete the integrals 
where data is lacking. Not all of these uncertainties apply 
to each of the numbers in the Table, and we indicate this by 
leaving blanks in the parenthesis. In the Table, the statisti- 
cal errors are always combined in quadrature; the systematic 
errors from g\ (the error in the low Q 2 bin is treated as due 
to uncertainty in gi, as it mostly is) are combined directly, 
and then combined in quadrature with those from the F 2 
terms; and the modeling errors in Ai and A2 are separately 
combined directly, and then combined in quadrature with 
each other. 

Further combining the statistical, systematic, and mod- 
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eling errors in quadrature gives the result 
A pol = 1.88 ±0.64 ppm. 



(21) 



For other choices of ^(Q 2 ), the resulting changes in A po i 
are small compared to the overall error limit quoted above, 
as may be seen in the A po i column of TableQU The current 
result differs from our previous 1.3±0.3 result [25j based on 
earlier data and a less sufficient treatment of the systematic 
errors. Other A po i results incompatible with zero are the 
2002 Faustov and Martynenko [27| value of 1.4 ± 0.6 ppm 
and the 2006 Faustov, Gorbacheva, and Martynenko [28j 
value of 2.2 ±0.8 ppm. 



C. Ordinary hydrogen Zemach and recoil corrections 

Form factor measurements have improved in the past 
decade largely due to the use of polarization transfer tech- 
niques [29] and to an understanding of how two- photon 
corrections impact the Rosenbluth measurements [30]. An- 
alytic form factor parameterizations new within the past year 
are available from Arrington, Melnitchouck, and Tjon 1 3lQ . 
who fit over all Q 2 where there is data, and from Arring- 
ton and Sick [Hj], who concentrate on the lower Q 2 data. 
The Zemach contributions from these, and two slightly older 
fits [H, [34j|, are listed in the third and fourth columns of 
TableQU One notices that modern form factors give larger 
radii and larger magnitude Zemach corrections than the old 
dipole form. 

(The fits of Ref. [32] are valid only for Q 2 < 1 GeV 2 , 
and for Q 2 above this value we supplement them with form 
factors taken from Ref. |3l|. The integrals are strongly 
weighted to lower Q 2 , so that if we supplement them with 
the dipole forms instead, the results would be the same to 
the number of figures given.) 

Using the same form factors, we quote the recoil correc- 
tions in the fifth column of TableQT] The bulk of the result 
comes from the one-loop corrections of Eq. (|15p . We also 
included a 0.46 ppm two-loop recoil correction from Bod- 
win and Yennie [15]. The latter are 0{a 2 ) beyond the Fermi 



TABLE II: Zemach radii, A z including 5^ d , and the recoil 
corrections, for four modern form factors. The dipole form 
factor is included only as a benchmark. The "target" As is 
—32.77 ± 0.01 ppm; the errors on As in the Table above are 
typically ±0.7 ppm. 



Form factor 


rp 


rz 


A z 


A?; 


Apoi 


As 




(fm) 


(fm) 


(ppm) 


(ppm) 


(ppm) 


(ppm) 


AMT [31J 


0.885 


1.080 


-41.43 


5.85 


1.88 


-33.70 


AS [32] 


0.879 


1.091 


-41.85 


5.87 


1.89 


-34.09 


Kelly [33] 


0.878 


1.069 


-40.99 


5.83 


1.89 


-33.27 


FW [34] 


0.808 


1.049 


-40.22 


5.86 


2.00 


-32.36 


dipole 


0.851 


1.025 


-39.29 


5.78 


1.94 


-31.60 



energy scale and are given by, 
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This correction is evaluated only in the static limit. For the 
0(a) recoil correction (leading order or LO), evaluating in 
the static limit gave a different sign and about a factor of 2 
smaller magnitude than using physical form factors. How- 
ever, this was possible only because the static evaluation of 
the LO is unexpectedly small due to striking internal can- 
cellations. To wit, the LO static correction evaluated with 
the measured k p is about 15 times smaller than it would be 
using K p = 0. Similar internal cancellations do not occur 
in the next-to-LO corrections, and we should not expect — 
albeit this is not verified — that evaluation with physical form 
factors would occasion big changes in the two-loop correc- 
tions. We also include an additional 0.09 ppm radiative 
correction noted by Karshenboim [l6| . 

An overall summary of calculated results for the ordinary 
hydrogen hfs is given in Table QTTJ along with, in the first 
line, the experimental value of the corrections in units of 
the Fermi energy. 

Regarding error limits for the Zemach correction, three of 
the four modern form factor fits give uncertainties in their fit 
parameters that allow an estimate of the uncertainty in Az 
obtained from the respective fits. There are, of course, cor- 
relations. For example, some of the data are cross sections, 
so that if the extracted Ge goes up, then Gm goes down. 
We estimated the uncertainty in Az by letting Ge vary to 
the maximum allowed by the respective authors's error lim- 
its, and doing so leads to variations in Az of ±0.085 ppm, 
±0.33 ppm, and ±0.80 ppm for the AS (H, Kelly |33|, 
and FW [34[ fits, respectively. This may be an argument 
for favoring the AS fit. However, the variations among the 
results for the different form factor fits are larger than some 
of the uncertainties just quoted, and we have taken the ap- 
proach of using the result from the AMT fit [3l| with an 
uncertainty chosen to accommodate the two most modern 
of the other fits. A similar choice has been made for the 
recoil corrections. 

The total of the hfs corrections gives a result that is 0.85 
ppm short of the data, with a quoted uncertainty of 0.78 
ppm. The goal of a 1 ppm calculation appears to have been 
reached, with the theory versus data difference barely over a 
standard deviation. There is no evidence of missing physics 
at this level. Also, the uncertainty in the polarizability term 
is now comparable to the uncertainty in the Zemach term, 
which is purely dependent upon the elastic form factors. 
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TABLE III: Summary of corrections for electronic hydrogen; 
Az, A^,, and A po i come from Tables HI and ITT1 



Quantity 


value (ppm) 


uncertainty (ppm) 


{E Ma (e-p)/E p F ) - 1 


1 103.48 


0.01 


Aqed 


1 136.19 


0.00 


A P I A P I A P 
"^vp ' "hvp ' "weak 


0.14 




A z (using [31]) 


-41.43 


0.44 


A p R (using [31]) 


5.85 


0.07 


Ap i (this work, using [3 11^ 


1.88 


0.64 


Total 


1102.63 


0.78 


Deficit 


0.85 


0.78 



D. Muonic hydrogen structure-dependent corrections 

For muonic hydrogen we, of course, keep me 7^ 0. There 
are no poles in the integrands of Eqs. ©I the numerators of 
the $i are zero when the denominators are. For numerical 
purposes, one can analytically divide to obtain 



-2t ( t 



vv+t 



- 1 



{2nr + 4/7*7) 
{2ut - At) , (23) 



Evaluating A* for the < Q 2 < Q\ = 0.0452 GeV 2 data 
gap is somewhat different from the electron case. Now, 
re = Tjx and is small (in the range to about 1) rather than 
very large, although t is still fairly large. A numerically 
good approximation for these ranges is 



(24) 



This leads to 
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where 70 is the forward spin polarizability, 

2a f°° dv 



To 



(25) 



(26) 



"V Jvth 



and v t h = TOtt + {fn 2 + Q 2 )/(2m p ). From data, 70 = 
[-1.01±0.08(stat)±0.10(syst)] x 10" 4 fm 4 [H, and c x = 
4.50 ± 0.35 ± 1.42 is from the analog of Eq. (jTH) but for 

ri(Q 2 )- J^dx 9l (x,Q 2 ) (H. 



TABLE IV: Contributions to A po i for muonic hydrogen. As 
in Table HJ statistical, systematic, and modeling errors are 
given in the parentheses. 



Term 


Q 2 (GeV 2 ) 


From 


Value w/AMT [31J F 2 


Ax 


[0, 0.0452] 


F2 and gi 


0.86(0.17)(0.67) ( ) 




[0.0452, 20] 


F 2 


6.77 ( ) (0.21) ( ) 






<h 


0.18(0.18)(1.62)(0.64) 




[20, 00] 


F> 


0.00 ( ) (0.00) ( ) 






9i 


0.11 () () (0.01) 


total Ai 






7.92(0.25)(2.30)(0.66) 


A> 


[0, 0.0452] 


92 


-0.12 () () (0.12) 




[0.0452, 20] 




-0.29 ( ) ( ) (0.29) 




[20, 00] 


92 


-0.00 ( ) ( ) (0.00) 


total A 2 






-0.41 () () (0.41) 


Ai + A 2 






7.51 (0.25) (2.30) (0.77) 


Apoi (ppm) 






351.( 12. )(107.)(36. ) 



The evaluation of the polarizability corrections for higher 
Q 2 is similar to the evaluation in the electronic case, and 
depends upon the same combination described previously 
of EG1 data for gi/Fi, Christy and NMC/E143 fits for F lt 
and supplements from Simula et al. and E155 fits where 
there is no EG1 data [H, 0, M. 112, HI ■ 

Table [jV] shows the breakdown of contributions to A po i 
for muonic hydrogen using the AMT [3j| elastic form factor. 
Error limits in the Table are combined the same way as for 
Table [D Finally combining the statistical, systematic, and 
model dependent errors in quadrature yields 



A pol = 351 ± 114 ppm. 



(27) 



Results for A po i using other form factors are, as in the elec- 
tron case, not greatly different on a scale set by the current 
systematic errors on A po i. Results are shown in Table [Vl 

Also in Table [V]are results for the Zemach and recoil cor- 
rections in the muon case. The structure dependent correc- 
tions become large compared to the electron case since they 
are, unlike the QED corrections, proportional to the lepton 
mass. The Zemach corrections follow simply from scaling 
the electron case with the new reduced mass. The recoil 
corrections are easily recalculated and include the two-loop 
corrections of Bodwin and Yennie quoted in Eq. (j22|) . The 
latter scale directly with the lepton mass; they were 0.46 
ppm for the electron case and are here 96 ppm. The extra 
radiative recoil corrections that accounted for 0.09 ppm in 
the electron case have been omitted. The vacuum polariza- 
tion part of these corrections are easy to scale to the muon 
case (l6j . but formulas are not available for the self energy 
part. These corrections were, for electronic hydrogen, small 
compared to the current overall accuracy of the final result. 
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E. Estimates regarding A2 

In this subsection we reconsider A2, first using the 
Wandzura-Wilczek approximation [36], and then consider- 
ing what existing proton data can tell us about the non- 
Wandzura-Wilczek part of A2. After reconsideration, we 
shall still believe that A 2 based on the EG1 model, with 
100% error limits, is satisfactory for the present hfs accu- 
racy goal. 

The structure function gi can be divided into 

g 2 (x, Q 2 ) = g¥ w (x, Q 2 ) + Q 2 ) , (28) 
where the Wandzura-Wilczek relation states 

gY W (x, Q 2 ) = -*(», Q 2 ) + r 9 -^^dy . (29) 

Jx V 

Although at high Q 2 , g^ w is the leading twist contribution 
and g 2 is the higher twist component, this formal division 
is well-defined all the way to Q 2 = 0. Hence, we can define 



B 2 (Q 2 ) = Bf w 



B 2 



BY W (Q 2 ) = / dxf3 2 (r)gY W (x,Q 2 ), 
Jo 



B 2 (Q 2 



dx0 2 {r)g 2 (x,< 



Substituting and manipulating yields 



dx0 i (j)g 1 {x,Q 2 ) 



(30) 



(31) 



in which 



&(t) - 4Vr(r + l)-4r-2Vfln (^^^) , (32) 

for the electron case. The function 03 has limits 

3 (t) = (Inr-21n2 + 4)VT-4r, r -> 0, 
0s(r) = -1- + ^ + ..., r^oo. (33) 



TABLE V: For muonic hydrogen hyperfine splitting: Zemach 
radii (as before, included for completeness), Az including 
5™ , recoil corrections, polarizability corrections, and the 
summed structure dependent corrections As, for four mod- 
ern form factors, with the dipole form factor included as a 
benchmark. Typical errors on As are ±120 ppm. 



Form factor 


rz 


A z 




Apoi 


As 




(fm) 


(ppm) 


(ppm) 


(ppm) 


(ppm) 


AMT [31J 


1.080 


-7703. 


931. 


351. 


-6421. 


AS [32] 


1.091 


-7782. 


931. 


353. 


-6498. 


Kelly [33J 


1.069 


-7622. 


931. 


353. 


-6338. 


FW [34J 


1.049 


-7479. 


939. 


370. 


-6170. 


dipole 


1.025 


-7311. 


935. 


362. 


-6014. 



The result is AJf w = -0.71 ±0.08 ±0.10± 0.01 using the 
same data and techniques as for determining B\ and Ai. 
For the muon case, one replaces 0a by 



03(r,n) = -0 2 (t, r e ) - 2r< + 2y/n{n + 1) 
+ 2^n(re + l) In (- 



2Vrln 



V7TT + 1 



(34) 



and obtains A^^p) = -0.57 ± 0.06 ± 0.10 ± 0.01. 

Notice that we obtain negative A 2 , which means that 
the main support in the integrals comes from regions where 
g 2 (x,Q 2 ) is positive. Hence, the main contribution to the 
hfs from the g 2 terms comes from the resonance region, and 
specifically from the region of the A(1232) resonance, since 
the existing data shows g 2 is positive there and negative 
elsewhere. As the integrals also have stronger support at 
low Q 2 , one specific need for more data would be in the 
higher resonance and continuum regions at low Q 2 . 

There is interest among hadronic physicists in the higher 
twist component g 2 (x,Q 2 ). In particular, the higher twist 
coefficient defined by 



d 2 {Q 2 



3 / x g 2 (x,Q )dx 
Jo 



x 2 {2g 1 + 3g 2 )dx 
(35) 

has been studied by experimenters. 

Osipenko et al. (3?} quote results for d 2 (Q 2 ) that are 
small at low and high Q 2 , but significant and positive for 
values within a decade on either side of Q 2 = 1 GeV 2 . 
However, they also give systematic errors, and these are 
very large. For comparison, Kao et al. [38] also model d 2 , 
and do not find a fixed sign. 

At one value of Q 2 , namely 1.3 GeV 2 , there is good data 
on the proton's g 2 (x,Q 2 ) from the RSS collaboration (39| . 
This data shows that g 2 is typically about half of g\ vw , 
at this Q 2 . If this be generally true then A2 is, of course, 
about half the A^ w values just quoted. 

We conclude that, partly because of the smallness of the 
contribution, the existing data allows us to use the A2 val- 
ues quoted in our Tables, with confidence that the generous 
percentage error limits will include any changes that will 
come with better data. Of course, we do want more com- 
plete data (there is already more complete data for neutron 
targets, extracted from polarized 3 He targets) and know 
that experimenters are also interested. 



III. CONCLUSIONS 

Our new result for the polarizability corrections to the hfs 
of the hydrogen ground state is, 



A pol = 1.88 ±0.64 ppm, 



(36) 



where the error limit includes both statistical and systematic 
uncertainties. The main new ingredient in this determina- 
tion is the recently published data on the spin-dependent 
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structure function giiv.Q 2 ) from the JLab EG1 collabora- 
tion [19(. The result is somewhat larger than our previously 
published A po i (25| , and the quoted uncertainty limit is also 
larger due to a better comprehension of the systematic error. 

There are also recent results on A po i from Faustov, Gor- 
bachova, and Martynenko [2£|. They have not used the 
EG1 data, relying instead on theoretically motivated fits to 
earlier data. Their result is somewhat larger than ours, but 
compatible within error limits. 

A consequence of the slightly larger A po i and larger un- 
certainty is that the calculated hfs is just within 1 ppm, and 
just about one standard deviation, from the experimental 
value. 

The goal of a 1 ppm hydrogen hfs calculation appears to 
be realized. One needs to make this claim with some diffi- 
dence; a well-known [15] paper is sometimes read as having 
made this claim in 1988. However, the claim only referred 
to the accuracy of the methods, and the authors themselves 
pointed out that the polarizability correction at that time 
was known only as compatible with zero to the 4 ppm level, 
and that the dipole form factor they used already differed 
systematically from the data. Indeed, Ref. [16] found that 
a better low Q 2 fit to the form factor data changed the 
Zemach contributions by about 2 ppm. Now, both form 
factors and structure functions are better known, and a 
claim of 1 ppm accuracy is plausible, with, at this level, 
no unknown terms remaining to be included. 

Further improvement in the calculation of hfs using elec- 
tron scattering data depends upon further improvement in 
the data and/or its analysis. The largest uncertainty cur- 
rently follows from systematic uncertainties in the inelastic 
structure functions. In Tables [T] and [jV] we give error limits 
separately for the statistical, systematic, and modeling er- 
rors. "Systematic" here means systematic errors only from 
the data, and "modeling" is the uncertainly estimated from 
the models that we use to complete the integrals where there 
is no data. The statistical errors are small. The largest er- 
rors are the systematic ones, which can only be improved 
by understanding the apparatus better, or by improved ap- 
paratus. 

In addition, the uncertainty in the Zemach term, which 
depends upon elastic form factors, is also noticeable. Even 
restricting to modern form factor fits, there is a 2% variation 
in the charge radius, with the fit [H, [4(| arguably most at- 
tentive to the low Q 2 data giving the largest result. Progress 
may depend not only on experimental progress, but also on 
a clearer understanding of the corrections needed to con- 
nect electron scattering cross section and polarization data 
to the form factors, and an assessment of how these cor- 
rections are implemented in current and future form factor 



fits. Atomic determinations, based upon Lamb shift mea- 
surements, are in line with the lower values and have about 
1% error limits [41]. Further, the atomic determinations of 
the charge radius will become remarkably more precise with 
the hoped for success of the muonic hydrogen Lamb shift 
experiment Q. 

The structure function giiy, Q 2 ) is less important for the 
hfs than g\, because the auxiliary function that multiplies 
it tends to be numerically small. This is good, given that 
cj2 is harder to measure than g-y, and there is little data for 
the proton. We included the 92 contributions with 100% 
uncertainty limits and we believe these suffice. The overall 
uncertainty still comes mainly from g\. We did give some 
further consideration to the 92 contributions, in particular 
calculating the Wandzura-Wilczek part and discussing the 
remainder. And certainly, more g2(x,Q 2 ) data would be 
welcome, to ensure that there are no surprises in, for exam- 
ple, the low Q 2 region at low x (where the higher resonances 
and continuum contribute). 

We have also given results for the structure dependent 
contributions to muonic hydrogen hfs. The total correc- 
tion is unambiguous, but for the massive lepton case the 
protocols for separating the recoil and polarizability terms 
seem not yet standardized. One hence needs to be watch- 
ful when adding together terms from different sources. We 
have quoted algebraically all the structure dependent terms 
to leading correction order (order ax mass ratio x Fermi 
energy) in the body of this paper, with matching conven- 
tions. Our result for the polarizability term is 

A P oi(M~p) = (351 ± 114) ppm , (37) 

and for the structure dependent terms overall, with the 
AMT elastic form factors, 

AsOu-p) = A Z +A P R + A po i = (-6421±140) ppm. (38) 

These are the biggest corrections for the muonic case. The 
QED corrections, in particular, are very nearly the same size 
as in the electronic case, since they do not have the mass 
proportionality that the structure dependent terms have. 
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